Systematics of low-lying states of even-even nuclei in the neutron-deficient lead region 

from a beyond-mean-field calculation 
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Background Nuclei located in the neutron-deficient Pb region have a complex structure, rapidly evolving as a function of 
neutron and proton numbers. The most famous example is 186 Pb where the three lowest levels are + states, the two 
excited + ones being located at low excitation energy around 600 keV. Coexisting structures with different properties 
are found in the neighboring nuclei. Many experiments have been performed over the last few years in which in-band 
and out-of-band 7 transition probabilities have been measured. 

Purpose A detailed interpretation of experimental data requires the use of a method going beyond a mean-field approach that 
permits to determine spectra and transition probabilities. Such methods have already been applied to selected isotopes 
in this mass region. Our aim is to provide a systematic investigation of this mass region in order to determine how well 
experimental data can be understood using a state-of-the-art method for nuclear structure. 

Method The starting point of our method is a set of mean-field wave functions generated with a constraint on the axial 
quadrupole moment and using a Skyrme energy density functional. Correlations beyond the mean field are introduced by 
projecting mean-field wave functions on angular-momentum and particle number and by mixing the symmetry restored 
wave functions as a function of the axial quadrupole moment. 

Results A detailed comparison with the available data is performed for energies, charge radii, spectroscopic quadrupole 
moments, E0 and E2 transition probabilities for the isotopic chains of neutron deficient Hg, Pb, Po and Rn. The 
connection between our results and the underlying mean field is also analyzed. 

Conclusions Qualitative agreement with the data is obtained although our results indicate that the actual energy density 
functionals have to be improved further to achieve a quantitative agreement. 

PACS numbers: 21.10.-k, 21.10.Re, 21.60.Jz, 27.70.+q, 27.80.-f-w 
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I. INTRODUCTION 

The study of the low-energy spectrum of neutron- 
deficient nuclei around 186 Pb has unveiled a rich variety 
of collective levels. In particular, several + states coex- 
ist at low excitation energy, which are in several cases the 
band head of a rotational sequence. This experimental 
result is interpreted as a manifestation of "shape coex- 
istence"; i.e., a condition where states at similar excita- 
tion energy have distinctly different intrinsic shapes [1—5]. 
Many data support this interpretation. Direct evidence 
is, indeed, provided by the systematics of a-decay fine 
structure [4, 6, 7], by the moments of inertia of the rota- 
tional bands built on several + states [4], by in-band and 
intraband transition probabilities [8], by g- factor mea- 
surements [9, 10] and by charge radii [11, 12]. 

Let us first recall how the low-energy spectrum of the 
Pb isotopes around N = 104 is interpreted in terms of 
deformation. For all isotopes, many experimental evi- 
dences indicate that the ground states are predominantly 
spherical [5] . The two lowest excited + states have been 
associated with two deformed structures, an oblate one 
lower in energy in the isotopes above 188 Pb, and a pro- 
late one that is lower in lighter Pb nuclei. The crossing 
of these two structures in the excitation spectrum leads 
to the unique situation where the three lowest levels in 
186 Pb are 0+ states [7]. 

By contrast, the ground states of the Hg isotopes down 



to i8U Hg arc interpreted as being oblate and weakly de- 
formed with /3 w -0.15 [13, 14]. From 198 Hg down to 
190 Hg, the yrast states for a given angular momentum 
have almost constant excitation energies and are inter- 
preted as the members of a rotational band based on 
the oblate ground state [15]. This simple pattern is dis- 
torted for the lighter even isotopes through the intrusion 
of a strongly deformed prolate band with (3 w 0.25 [3, 4]. 
The excitation energies of the states in the prolate band 
evolve in a nearly parabolic manner as a function of the 
neutron number, cf., for example, Fig. 1 in Ref. [16] or 
Fig. 10 in Ref [5]. 

The production rates of light Po isotopes are much 
smaller than those of the corresponding Hg and Pb iso- 
tones; hence, less experimental data have been collected, 
but they give a clear indication for an even more com- 
plex evolution of their structure as a function of neu- 
tron number. The heavier Po isotopes down to 196 Po 
have near-spherical ground states [5]. For 188 - 192 p , the 
analysis of their a-decay fine structure indicates that the 
ground-state wave function contains a significant contri- 
bution from deformed configurations [17]. Similar con- 
clusions have been drawn from charge radii that indicate 
an increasing softness of the nucleus against quadrupole 
deformations with decreasing neutron number [18]. 

Our study will be limited to even-even nuclei. Let us, 
however, mention that data for odd-mass nuclei, where 
the low-lying spectrum can often be interpreted by the 
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coupling of one particle to the low-lying states of the ad- 
jacent even-even nucleus, corroborate the interpretation 
of the spectral data obtained for the light even-even iso- 
topes of Hg, Pb, Po and Rn [5]. 

There are two different approaches to interpret the 
complex and rapidly changing structure of nuclei in this 
mass region. One is based on the shell model, where 
the emergence of low-energy intruder states is interpreted 
as resulting from proton excitations across the Z = 82 
closed shell [3, 19, 20]. Although the shell model in- 
cluding np-mh excitations across the shell gaps pro- 
vides a consistent and microscopic description of the phe- 
nomenon of shape coexistence throughout the nuclear 
chart [5], the high dimensionality of the corresponding 
model space renders calculations for heavy open-shell nu- 
clei untractable. For such nuclei, however, the np-mh ex- 
citations can be handled within the algebraic framework 
of the interacting boson model, which truncates the shell- 
model space through the approximation of J = and 
J = 2 coupled nucleon pairs treated as bosons [21, 22]. 

An alternative and intuitive approach is based on 
mean-field methods, which provide a description of these 
structures in terms of shapes and deformed shells. In this 
framework, the low-lying spectrum is explained by the 
presence of multiple local minima in the deformation en- 
ergy surface. These minima can in turn be related to gaps 
in the single-particle spectrum. Although shape coexis- 
tence can be interpreted within pure mean-field models, 
a detailed description of the spectrum requires to "go be- 
yond the mean-field" by means of symmetry restoration 
and taking into account fluctuations in the deformation 
degrees of freedom. The beyond self-consistent mean- 
field method that has been developed by several groups 
combines the projection techniques with the generator 
coordinate method (GCM) [23-29]. 

Throughout this work, we use the axial quadrupolc 
moment as a collective variable. Mean-field wave func- 
tions are constructed by Hartrec-Fock (HF)+BCS calcu- 
lations. They are projected on angular-momentum and 
particle number to form a basis of states that are mixed 
by the GCM. The interaction used here is the SLy6 para- 
metrization of the Skyrme interaction that describes well 
the systematics of low-lying states of neutron-deficient 
Pb isotopes [30, 31]. Selected results for some of the nu- 
clei discussed here have already been published earlier in 
Refs. [14, 18, 32-34]. 

The paper is organized as follows. In Sec. II, we will 
give a brief outline of the method used to calculate the 
spectroscopic properties of low-lying states. In Sec. Ill, 
the calculated deformation energy curves, low-energy ex- 
citation spectra, charge radii, kinetic moments of in- 
ertia, and electric monopole and quadrupole transition 
strengths are presented and discussed in comparison with 
results of our previous study of Pb isotopes and the avail- 
able data. A summary of our findings and conclusions is 
given in Sec. IV. 



II. THE METHOD 

As in our previous works [31], the starting point 
of our method is a set of HF+BCS wave functions 
\q) generated by self-consistent mean-field calculations 
with a constraint on the axial mass quadrupolc moment 
q = (q\2z 2 — x 2 — y 2 \q) . Dynamic correlations associated 
with symmetry restorations and fluctuations in the shape 
degree of freedom are introduced by particle-number and 
angular-momentum projection in the framework of the 
GCM. Limiting ourselves to axially symmetric configu- 
rations, the final wave function for the correlated state 
\JM;k) is given by the superposition of symmetry re- 
stored mean-field wave functions 

\JM;k) = J2fk(q)\JMq), (1) 

9 

where k ~ 1, 2, . . . labels different collective states for a 
given angular momentum J. The variable q is the generic 
notation for the deformation parameters. The symmetry 
restored mean-field wave function is constructed as 

\JMq) = ^—J2g-kP{iKP z P N \q), (2) 

■Nj,M,q x 

where Mj.M.q is a normalization factor. Pmk projects 
out eigenstates of J 2 and J z with eigenvalues H 2 J(J+ 1) 
and KM in the laboratory frame or HK in the intrin- 
sic frame, respectively, whereas P N and P z project out 
eigenstates of the particle-number operator for neutrons 
and protons with eigenvalues N T , t = n, p. 

The weight functions and the energies of 

the states \JM;k) are the solutions of the Hill- Wheeler- 
Griffin (HWG) equation [35] 

J2 [n J (q, q') - E( M J (q, q')] //(</) = . (3) 
q' 

The ingredients of Eq. (3) are the norm kernel 
N J (q,q') = (JMq\JMq') and the energy kernel 
U J (q,q'), which in our calculation is given by a multi- 
reference energy density functional that depends on the 
mixed density matrix [36] . The formulae used to evaluate 
the energy and the norm kernels have been presented in 
Rcf. [24]. 

The weight functions f£ (q) in Eq. (1) are not orthogo- 
nal. A set of ortho normal collective wave functions gf. (q) 
can be constructed as [37, 38] 

9Mq)=J2(^ J ) 1/2 ^l')fk(9')- (4) 
9' 

It has to be stressed, however, that the \g.jk(q)\ 2 quantity 
does not represent the probability to find the deforma- 
tion q in the GCM state \JM;k). In addition, in the 
absence of a metric in the definition of the correlated 
state \ JM;k), Eq. (1), the values of gjk(q) for a con- 
verged GCM solution still depend on the discretization 
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chosen for the collective variable q, which is not the case 
for observables such as the energies or transition proba- 
bilities. 

Since the correlated states \JM;k) have good angu- 
lar momentum, their spectroscopic quadrupole moments 
Q s (Jk) 



Neutron 



Qs(Jk) = o —j 

xE/fW W\\Q2\\Jq)f^(q) (5) 
g'.g 

and the reduced electric quadrupole (E2) transition 
strengths between them 

B{E2-J ki -> J kf ) 

q^lE^/V) WII&IIW£(«)| 2 (6) 



2J, 



9', 9 



are calculated directly in the laboratory frame without 
approximation. The reduced matrix elements entering 
both expressions are determined as 



(Jfq'\\Q2\\Jiq) 

(2J / + l)(2J i + l) 



-2 

E 

Af=-2 



M —M 



Jo 

(7) 

where Q2M = er 2 Y2M is the electric quadrupole moment 
operator. To relate the moments in the laboratory frame 
to intrinsic deformation parameters, one can define two 
dimensionlcss quadrupole deformations in the same way 
as in the static rotor model. The first one, (3^ , is related 
to the spectroscopic quadrupole moment (for states in a 
K = band) by: 
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2J + 3 
J 



Qs(Jk) (8) 



and the second, , is related to the reduced E2 transi- 
tion strength 



P (t) {J h ki -> Jf,k f ) 



4tt B(E2;Ji,k t ->J f ,kf) 



3ZR 2 



e 2 {Ji020\JfO)< 



(9) 

The radius R appearing in both expressions is given by 
R = 1.2A 1/3 fm, with A being the mass number and 
(<7i020|J/0) a Clebsch-Gordan coefficient. The nuclear 
matrix element entering the electric monopole decay from 
\JM;k) to \JM\k') through the emission of conversion 
electrons is determined by 



{JM;k'\eJ2 p r 2 p \JM ] k) 



5R 2 



(10) 




FIG. 1: (Color online) Nilsson diagram of the eigenvalues 
of the single-particle Hamiltonian for neutrons in 190 Po as a 
function of the axial deformation parameter /3, obtained with 
the Skyrme interaction SLy6, see Eq. (11). 



where p is an index running over all proton single-particle 
states. Since the electric transition matrix elements are 
calculated in the full model space of occupied single- 
particle states, there is no need to introduce effective 
charges, and we use the bare charge for protons instead. 

The parametrization SLy6 [39] of the Skyrme inter- 
action and a density-dependent zero-range pairing force 
with the same strength of —1250 MeV fm 3 for neutrons 
and protons and a soft cutoff at 5 MeV above and be- 
low the Fermi energy as defined in Ref . [40] , are adopted 
in the construction of mean-field wave functions and the 
configuration mixing calculations. As required by the 
SLy6 parametrization, the full two-body center-of-mass 
correction is included in the variational equations to gen- 
erate the mean field and in the calculation of the pro- 
jected GCM energies. 



III. RESULTS AND DISCUSSION 

A. General comments 

In mean-field models, self-consistent and non-self- 
consistent ones alike, there is an intimate relation be- 
tween the total binding energy and the density of levels 
around the Fermi energy in the Nilsson diagram of single- 
particle energies. This is obvious for non-self-consistent 
macroscopic-microscopic models, where the extra bind- 
ing from shell structure is explicitly calculated as a "shell 
correction" to the total binding energy, and, therefore, 
can be easily isolated, but the same mechanism is hid- 
den in the total energy of self-consistent mean-field mod- 
els. A lower-than-averagc density of single-particle levels 
around the Fermi energy results in extra binding, whereas 
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a larger-than-average value reduces binding. Therefore, 
large gaps near the Fermi energy in the Nilsson diagram 
often correspond to minima in the deformation energy 
curve, whereas a large bunching of levels close to the 
Fermi energy usually corresponds to barriers between 
such minima. The Fermi energy does not have to be 
located exactly at the gap, but might be slightly below 
or above it. With increasing distance of the Fermi energy 
from the gap, however, the shell effect can be expected 
to become less pronounced. In any event, shape coexis- 
tence is the fingerprint of a variation of shell structure 
around the Fermi energy that opens up or closes gaps 
with deformation. 

To have near-degenerate collective states with the same 
quantum numbers, but different deformations, there has 
to be a mechanism that prevents their mixing. As states 
of different deformations arc a priori not orthogonal, the 
diagonalization of the Hill- Wheeler- Griffin equation (3) 
can easily produce collective wave functions spread over 
a large range of deformations. Large non-diagonal en- 
ergy kernels will amplify this effect. Therefore, the sole 
presence of a large barrier between the minima is not suf- 
ficient to prevent mixing. To achieve a weak coupling be- 
tween mean-field configurations, the matrix elements of 
the overlap and the energy kernels between them must 
be small. A mechanism leading to such suppression of 
the matrix elements is provided by intruder states; i.e. 
single-particle levels that are downsloping with deforma- 
tion from the next major shell and, therefore, have a 
parity opposite to that of the levels at the Fermi en- 
ergy Assuming conserved parity and neglecting pairing 
correlations, two HF states with a different number of 
positive-parity and negative parity single-particle levels 
do not overlap and are weakly coupled in the GCM. This 
is no longer the case when pairing correlations are taken 
into account, but one can expect that the overlap between 
the mean-field wave functions remains small. 

The detailed results of this extensive study are avail- 
able in the supplemental files tabl.dat and tab2.dat. 
They include excitation energies, charge radii, spec- 
troscopic quadrupole moments, electric monopole and 
quadrupole transition strengths for the low-lying states 
that were computed for total angular momentum val- 
ues up to J = 10ft, in the neutron-deficient 17 ^- 1( ^ A Yig. 

180-194p b; 186-210p o and 194-204 Rn isotopes . We focus 

here on the evolution of low-lying states with neutron 
and proton number and discuss the similarities and dif- 
ferences between these four isotopic chains. 



B. Nilsson diagrams 

Figures 1 and 2 display the Nilsson diagram of neutron 
and proton single-particle energies as a function of a di- 
mensionless quadrupole deformation parameter j3. This 
parameter is related to the intrinsic mass quadrupole mo- 
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FIG. 2: (Color online) Same as Fig. 1, but for the protons. 



ment q of the HF+BCS states by the relation 
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There is a large number of gaps visible in the Nilsson 
diagram for neutrons. First, there is a spherical one 
100, below the 1*13/2+ intruder level. There 
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are many deformed gaps of slightly varying size on the 
oblate side for every even neutron number from N = 102 
to N = 110 at increasingly large deformations up to 
(3 ~ —0.2. These gaps are all located between the spread- 
ing magnetic substates of the H13/2+ shell. There are also 
gaps of varying size for every even neutron number be- 
tween N — 98 and N = 108 on the prolate side. Due 
to the downsloping K = 1/2 level from the 2<7 9 / 2 + shell 
above and the upsloping levels from the lft.9/2- shell be- 
low, however, they correspond to configurations with a 
different number of occupied intruder levels. Because of 
several level crossings, the deformations of the gaps on 
the prolate side arc not ordered according to N. The 
largest ones are N = 106 for f3 ~ 0.3 and N — 104 for 
j3 ~ 0.5, respectively. 

For the protons, there is a large spherical shell gap at 
Z = 82, a gap at Z = 80 extending from sphericity to 
an oblate shape with (3 values down to —0.2 as well as 
a smaller oblate Z = 82 gap at j3 ~ —0.25. At normal 
deformation, there are no significant proton gaps on the 
prolate side that come close to the Fermi energy, except 
for a small Z = 86 one at /3 ~ 0.5. 

At the deformation corresponding to the deformed 
neutron gaps, the downsloping high-j proton intruder or- 
bitals have dropped to energies close to the Fermi energy, 
in most cases even well below, and, thus, have a sizable 
occupation. For this reason, some single-particle orbitals 
composing the corresponding deformed mean-field wave 
function are very different from those of the spherical con- 
figuration and the overlap between both is small. As dis- 
cussed above, a small overlap creates favorable conditions 
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FIG. 3: (Color online) Energy curves for the particle- 
number-projected HF+BCS states (N&cZ) for 176 - 194 Hg, 

180-194p b) 186-204p O5 194-204^ ag & function Q f ^ in _ 

trinsic quadrupole deformation parameter ft. All energies are 
normalized to the spherical shape (/3 = 0), but with an ad- 
ditional energy shift of 1 MeV between two neighboring iso- 
topes. The lowest configuration of a given nucleus in the 
deformation region of interest is indicated with a bullet. 



for a decoupling between configurations corresponding to 
different shapes. 



We have checked that the Nilsson diagrams of both 
protons and neutrons change only marginally for the nu- 
clei covered in the present study. 



C. Systematics of deformation energy curves 

The energy curves obtained after projection onto par- 
ticle numbers are presented in Fig. 3. The mean- field 
ground states for all neutron-deficient Pb isotopes are 
spherical. The same result has been obtained in cal- 
culations using the parametrization D1S of the Gogny 
force [41]. However, it can be modified by a lower pair- 
ing strength [42] or the use of a different pairing func- 
tional [43]. For all Pb isotopes, except 180_182 Pb, the 
curves also display an oblate minimum with (3 ~ —0.20. 
The excitation energy of this minimum decreases from 
0.76 MeV to 0.50 MeV for N going from 112 to 108, and 
increases again up to 1.63 MeV when decreasing the neu- 
tron number further down to 102. The well depth follows 
a similar evolution. It starts from 0.46 MeV for N = 112, 
is maximal at 0.93 MeV for N = 108 and decreases again 
down to 0.49 MeV for N = 102. On the prolate side, 
the energy curves present an inflexion point at j3 ~ 0.30 
for N > 110 and a minimum for lower N values down 
to 100. The excitation energy of this minimum decreases 
from 1.47 MeV in 190 Pb down to 0.61 MeV in 184 Pb and 
raises up again to 0.82 MeV in 182 Pb. 

The deformation of the prolate minima of the Pb iso- 
topes follows closely the neutron gaps in the Nilsson di- 
agram. Their deformation is fa — 0.32 for N = 100 in- 
creases slightly at N = 102, and then decreases up to 
N = 106. For 108 < N < 112, this prolate minimum 
disappears as the Fermi energy approaches a region with 
densely bunched levels in the Nilsson diagram. 

There are no spherical minima for the Hg, Po and Rn 
isotopes discussed here. The minima of the deformation 
energy curves for these nuclei arc not always in a one-to- 
one correspondence with the gaps in the Nilsson diagram. 
This is because the gap for one nucleon species is often lo- 
cated at a deformation corresponding to densely bunched 
levels for the other. Hence, the minima correspond to the 
most favorable compromise between the shell effects for 
neutrons and protons. 

The ground states that we obtain for the Hg isotopes 
are oblate above N = 106 and prolate below this neutron 
number. At N — 106, the oblate and prolate minima 
are nearly degenerate. The energy curves are very shal- 
low for N = 98 and 96. The deformation of the oblate 
minimum corresponds to the Z = 80 proton gap for all 
values of N, whereas the prolate minima of the isotopes 
with N > 98 correspond to the many gaps obtained in 
the neutron Nilsson diagram around j3 ~ 0.3. For the 
soft N = 96 isotope, the prolate minimum is shifted to 
a weakly deformed configuration with /3 ~ 0.13, close in 
energy with the weakly deformed oblate minimum. The 
oblate minimum is very rigid for N = 114 and becomes 
softer with decreasing N. Such a behavior can be related 
to the neutron single-particle levels close to the Fermi en- 
ergy. For N > 108, the density of neutron levels is high 
and their energy is rapidly increasing when (3 becomes 
more negative. The situation changes at N = 106 and 
below, where the density of neutron levels is low and the 
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FIG. 4: (Color online) Quadrupole deformation energy sur- 
face for 184 Hg, normalized to the absolute minimum and pro- 
jected on particle numbers. Each contour line is separated 
by 0.2 MeV. The inset shows the energy as a function of 7 
deformation along the path joining the two axial minima. 



energy varies more slowly with deformation. 

The energy curves of the Po isotopes evolve in a way 
similar to those of the Hg isotopes. In particular, the 
transition between oblate and prolate ground states ap- 
pears at the same neutron number. However, there are 
differences in the details. The oblate minimum is al- 
ways shallower, and the oblate deformation of the ground 
state increases steadily with decreasing N. For TV = 108, 
which corresponds to the lightest Po isotope with an 
oblate ground state, the deformation parameter /3 is sig- 
nificantly larger for Po than for Hg. By contrast, for 
102 < N < 106 the prolate ground states have a smaller 
deformation for Po than for the corresponding Hg iso- 
topes. 

The energy curves of the Rn isotopes evolve in a way 
similar to those of the Po isotopes, but with deeper oblate 
and prolate minima. 

Our calculations are limited to axial deformations. 
One cannot exclude, however, that in some cases the min- 
ima appearing in the axial deformation energy curves ac- 
tually correspond to saddle points of the energy surface in 
the full /3~7 plane. Several scenarios are possible. First, 
the axial minima can be separated by a barrier with the 
energy rising all the way from both sides, such that they 
are true minima. In Ref. [31], we have checked that this 
is the case for the spherical, prolate and oblate minima 
found in the 182 - 194 p D isotopes studied there. But there 
could be cases where there is no barrier between the ax- 
ial minima. A second possibility is that the energy rises 
monotonously from 7 = 0° to 60°. Then, the higher- 
lying minimum is in fact a saddle point. This happens 
for many well-deformed nuclei, and also for many 7-soft 
ones. A third case occurs when there is a triaxial mini- 
mum between the two axial ones. In this case, depending 
on the appearance of barriers, one or even both minima 
in the axial energy curve might be saddle points. 




FIG. 5: (Color online) Same as Fig. 4 for 200 Po. 

The Gogny interaction Dls very often gives deforma- 
tion energy surfaces very similar to those obtained with 
the Skyrme interaction SLy6 used here. A systematic 
survey of energy surfaces in the plane using Dls 
can be found in [44]. For the nuclei studied here, the 
study of [44] indicates that minima found in a calcula- 
tion restricted to axial shapes have to be suspected to be 
saddle points when the ground state and the secondary 
minimum in the axial energy curve are nearly degenerate 
and are obtained for similar values for \/3\. 

Two examples of particle-number projected energy sur- 
faces in the full /3—y plane obtained with SLy6 are dis- 
played in Fig. 4 for 184 Hg and Fig. 5 for 200 Po, respec- 
tively. The variation of the energy as a function of 7 for 
a path going from the oblate to the prolate minimum is 
shown as an inset in the figures. 

Fig. 4 illustrates the case where there is small bar- 
rier between the two minima. For 200 Po, the absolute 
minimum is obtained for a non-axial configuration, with 
f3 k, 0.1 and 7 « 10°. This nucleus is a rare case where 
both minima found in an axial energy curve actually cor- 
respond to saddles in the full (3-j plane [45] . The prolate 
saddle is at almost the same \(3\ value as the nearby tri- 
axial minimum and is located just about 100 keV higher, 
whereas the oblate saddle is at somewhat smaller \/3\ val- 
ues and higher excitation energy. The heaviest Po and 
Rn studied here can be expected to have deformation en- 
ergy surfaces of similar 7-soft topography, although not 
necessarily with a triaxial minimum. 



D. Effects of angular momentum projection on 
energy curves 

The energy curves projected on J = are given in 
Fig. 6. The energies are drawn at the intrinsic quadrupole 
moment (or the /?- value) of the mean-field wave function 
that is projected. This is the most convenient way to 
plot the results obtained after projection on angular mo- 
mentum. However, it has to be kept in mind that, after 
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FIG. 6: (Color online) Same as Fig. 3, but with projections 
onto both particle numbers and angular momentum J = 0. 

angular-momentum projection, it is only at large defor- 
mation and for J > that the (3 value of Eq. (11) provides 
a rough estimate for the values obtained through Eqs. (8) 
or (9) from observable quadrupole moments. At small 
deformations, the intrinsic deformation does not have a 
relation to an observable. Therefore, the interpretation 
of these energy curves requires some caution [31, 41]. 
The energy A_E rot = Ej =0 ((3) — E({3) gained from the 
restoration of rotational symmetry is displayed on the 
bottom panel of Fig. 7 for 186 Hg, 188 Pb, and 190 Po. The 
spherical configuration is purely a J = configuration; 
hence, projection on J = does not bring any gain of 
energy. A slight deformation of the mean-field is suffi- 
cient to introduce higher J components (or, in a shell- 
model language, particle-hole excitations) in the mean- 
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FIG. 7: (Color online) Effect of projection on J = on the 
topography of deformation energy curves at small deforma- 
tion. Top panel: energy curves after projection on parti- 
cle numbers (N&zZ) and on angular momentum J = as 
a function of /3 for 1S6 Hg, 18S Pb, and 190 Po. Lower panel: 
The corresponding energy gained from projection on J = 0; 
AE TOt (P) = E J=0 (l3)-E(l3). 



field wave function. Projection on J = gives then an en- 
ergy gain, which at small deformation is almost symmet- 
rical around the spherical point and increases rapidly (in 
absolute value) to reach about 3 MeV around \(3\ w 0.1. 
At larger deformation, the energy gain still increases fur- 
ther, but at a slower rate, as illustrated in Fig. 7. For 
nuclei with deformed minima, this generic behavior usu- 
ally makes the energy curve more rigid at small deforma- 
tion values and softens it at larger ones. For nuclei with 
a spherical minimum in the mean-field energy surface, 
such as the Pb isotopes, the minimum is shifted to small 
absolute values of the deformation (3 for all values of 7, 
forming a kind of "Mexican hat" around the spherical 
point. 

Neither the mean-field wave functions nor the wave 
functions obtained for a given set of quantum numbers 
after symmetry restoration are orthogonal. The overlap 
between J = wave functions is large for all small de- 
formation values of the mean-field. It turns out that the 
states "at the bottom of the mexican hat" are not only 
almost degenerate, but also have very large overlaps close 
to one, meaning that after projection they describe the 
same physical state. We will call a state in this very 
particular minimum a " correlated spherical state" . This 
state has the property to appear only for J = values, 
but not for higher angular momenta, and therefore, it is 
not the head of a rotational band. This allows to dis- 
tinguish the correlated spherical states from "truly de- 
formed" minima at small deformation that appear at all 
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J and lead to the appearance of a rotational band. In 
any event, at small deformation, the labeling of projected 
wave functions by a mean-field deformation has a very 
limited meaning. 

The main findings for the effect of projection on J = 
on the energy curves of the nuclei discussed here are: 

• Using the definition of a correlated spherical state, 
the ground state is spherical for all Pb isotopes. 
A prolate minimum around /3 ~ 0.3 is obtained for 
the lightest isotopes up to TV = 106. The oblate 
minimum that was obtained before projection on 
angular momentum for isotopes above N = 106 is 
not clearly visible anymore, as it merges with the 
spherical minimum. 

• For the Hg, Po and Rn isotopes, angular momen- 
tum projection does not modify the topography of 
the energy curves significantly. Both the prolate 
and the oblate minima become more bound with 
respect to the energy of the spherical point. For 
some nuclei, the minimum is also shifted to slightly 
larger deformation. 

• For Po isotopes with 112 < N < 120, the abso- 
lute minimum of the energy curves moves from an 
oblate to a prolate shape after projection on J = 0. 
These Po nuclei are the only ones where projection 
results in a change of the sign of the deformation. 

When triaxial shapes are included, angular momentum 
projection shifts the deformed minima from axial shapes 
to slightly non-axial ones [26-28], in the same manner 
as the spherical minimum is shifted to slightly deformed 
intrinsic shapes. For light nuclei, symmetry-restored 
GCM calculations including triaxial shapes have been 
performed and this shift does not qualitatively change 
the interpretation of the minima. Because of their high 
computational cost, such calculations have not been per- 
formed yet for nuclei as heavy as the ones discussed here 
within an energy density functional framework. 

While the mean-field ground states are at normal de- 
formation for all nuclei considered here, we have to note 
that after angular momentum restoration, the ground 
states of 186 > 188 Po and 194,196 Rn correspond to the projec- 
tion of a superdeformed configuration beyond the range 
of deformations displayed in Fig. 6. A similar result for 
the lightest Po and Rn isotopes is also obtained when 
using the Gogny Dls force and adding a rotational cor- 
rection [44]. This finding is the consequence of the en- 
ergy gain from angular momentum projection growing 
with deformation, cf. Fig. 7. It is in contradiction with 
the data and is an artefact of the low surface energy co- 
efficient a sur f = 17.7 MeV of the parametrization SLy6. 
We have checked that with the SLy4 parametrization, 
which has a slightly larger surface energy coefficient of 
a surf = 18-4 MeV, the deformation energy curves are 
stiffcr such that the superdeformed minima remain above 
the normal deformed ones after projection. On the other 



hand, at the mean-field level, the relative excitation en- 
ergies of the various minima in the Pb region are much 
better described by SLy6, as, overall, SLy4 gives excita- 
tion energies that are too large. A similar result has been 
found for 240 Pu in Ref. [29]. This points to the need to fit 
dedicated paramctrizations for bcyond-mean-ficld calcu- 
lations in the future. We have limited our configuration 
mixing calculations to deformations up to j3 « 0.6. We 
have checked that the normal- and superdeformed states 
are sufficiently decoupled that the low-lying states at nor- 
mal deformation discussed here are not affected by using 
these restricted configurations. 



E. Systematics of low-lying states 

1. Collective wave functions 

The collective wave functions g^.(q), as defined by 
Eq. (4), are spread over a large range of deformed mean- 
field wave functions. However, it is often still possible 
to classify them as spherical, oblate or prolate by look- 
ing at the dominant configurations. Such a classification 
is greatly helped by studying the mean deformation of 
the mean-field components of a collective wave function 
defined by 

P Jk ^Y,^)\9Jk{<l)\ 2 - (12) 

9 

This quantity has to be taken with a grain of salt. First, 
because of the non-orthogonality of the basis of projected 
mean- field states, \gjk{q)\ 2 does not represent the proba- 
bility to find a given mean-field state \q) in the projected 
GCM state \ JM;k). Second, /3jfc can not be expected 
in general to correspond to a deformation deduced from 
electromagnetic transitions or spectroscopic moments. 
Still, this quantity turns out to be useful when analyzing 
the mixing of oblate, spherical and prolate wave func- 
tions. For simplicity, we will label the collective states as 
oblate, prolate or spherical when the mean deformation 
/3jfc has a value close to the deformation of the mean-field 
state with the largest weight in the collective wave func- 
tion. When this is not the case, the state will be denoted 
as " spread" . 

The collective wave functions gjk(q), Eq. (4), of the 
lowest J = states are plotted in Figs. 8-11 as a function 
of the deformation /3 for i80-i94p bj i76-i94 Hg i86-204p o 

and 194 ~ 204 Rn, respectively. The third + is displayed 
only when it presents a special interest. The /?ofc-value 
is indicated for each state by the symbols »(0f ), ■(0^), 
and A(O^), respectively. 

For all Pb isotopes, the GCM ground state is predom- 
inantly spherical, in agreement with the data. The two 
lowest excited states are dominated by either oblate or 
prolate components. As we shall see below, the prolate 
state is lower in energy for the lighter isotopes up to 
N = 106 and the oblate one is lower above this value 
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FIG. 8: Collective wave functions for the first two and some 
selected third 0+ states in iso-^Pb. The 

mean deformations 

/3jk = ^2 q P(q) (<?0fc(<2)| 2 are indicated with the symbols »(0f ), 
■ (0+), and A(0+). 



U) 



U) 



0.4 
0.0 
-0.4 

0.4 
0.0 
-0.4 

0.4 
0.0 
-0.4 

0.4 
0.0 
-0.4 

0.4 
0.0 
-0.4 




=102H 

1 1 1 1 1 1 1 n 1 1 1 



N=106-^ 

I ' J ' | ' I -. ' I i I 





=104- 

1 1 1 1 1 1 1 .n 1 1 1 



N=108" 

I ' .I ■ I 




-0.4 0.0 0.4 0.8-0.4 0.0 0.4 0.8 

P P 

FIG. 10: Same as Fig. 8, but for 1S6 - 204 Po. 



_5 o.o 




U) 



0.4 
0.0 
-0.4 

0.4 
0.0 
-0.4 

0.4 
0.0 
-0.4 

0.4 
0.0 
-0.4 



N=96 

I i I i I i I i I i I i 




=iocq 

■ 1 1 1 1 1 n 1 1 1 



=104- 

K++-H-H-K-H 




N=98 

i I i I i I i I i I 



102i 

t i i . i i i i i n i i i 



N=106- 






-0.4 0.0 0.4 0.8-0.4 0.0 0.4 0.8 

P P 



FIG. 9: Same as Fig. 8, but for 176 " 194 Hg. 



while they are nearly degenerate for N = 106. Although 
they are not separated by a sizable barrier, the spherical 
and oblate states are well separated. The occupation of 
intruder single-particle states in the oblate configuration 
seems to prevent a mixing. Indeed, at /3 ~ —0.15, the 
downsloping K = 7/2 intruder orbital from the 2/ 7 / 2 - 
shell above the spherical Z — 82 shell closure crosses the 
upsloping 3s!/ 2 + orbital. In the same way, the prolate 
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FIG. 11: Same as Fig. 8, but for 194 ~ 204 Rn. 



state does not mix significantly with the spherical and 
oblate ones as its dominant contributions have a large 
occupation of several downsloping intruder levels on the 
prolate side. 

The mean deformation /3oii Eq. (12), for the ground 
states of 176 Hg and 178 Hg displayed in Fig. 9 is close to 
zero. Indeed, the small value of /3oi is the result of a 
cancelation between the contributions of oblate and pro- 
late configurations. This large mixing of different shapes 
can be related to flat energy curves given for these two 
isotopes in Fig. 3. For both isotopes, the first excited + 
state is prolate. 

In fact, from the discussion of triaxial energy sur- 
faces above, one can assume that both of these nuclei 
are soft in the 7 degree of freedom that is not explic- 
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itly considered here. However, we recall that in sym- 
metry restored GCM, prolate and oblate configurations 
are directly mixed by the Hill- Wheeler- Griffin equation, 
Eq. (3), independently of the height of a barrier on the 
axial path between them. Therefore, in an axial GCM 
calculation, as done here, one cannot distinguish the wave 
function of a 7-soft nucleus from that resulting from large 
mixing of two minima separated by a barrier. 

For the Hg isotopes 100 < N < 104, the prolate con- 
figuration is the ground state in contradiction with the 
data, as we shall see below, and the first excited + state 
is a spread state. The balance between oblate and pro- 
late configurations in the ground state is reversed above 
N = 104, with a larger weight on the oblate side for the 
heaviest isotopes represented in the Figure. The first ex- 
cited + state is predominantly prolate for N = 96, 98, 
106 and 108 and has a large spreading above 108. The 
second excited state is oblate for N = 106 and 108. 

For the Po isotopes, the ground states for N = 102 
and 104 are the only ones to be dominated by prolate 
configurations. The ground states of the heavier iso- 
topes have a large spreading. Figure 5 indicates that 
for these isotopes, the almost equal mixing of prolate 
and oblate configurations results from the 7 softness of 
their energy surface. The mean deformation [3jk of ex- 
cited states for 102 < N < 106 is small, reflecting the 
nearly equal weight of oblate and prolate configurations. 
From N = 108 to the heaviest Po isotopes studied here, 
the first excited + state is dominated by oblate com- 
ponents, whereas the second excited one is prolate for 
N = 108 and 110. Experimentally, the ground states of 
the Po isotopes are interpreted as being spherical down 
to N = 112. The lighter isotopes are inferred to posses 
deformed ground states [5]. 

The collective wave functions of the Rn isotopes behave 
in a way similar to those of their Po isotones. 
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FIG. 12: (Color online) Systematics of the excitation energy 
of the first three + states in the Hg, Pb, Po and Rn isotopes. 
Lines connect the states with similar dominant configurations. 
The available experimental data [16, 46, 47] for the low-lying 
+ states are plotted for comparison. 



isotopes, it has a very spread character. Also, 
for 100 < N < 104, the agreement between calcula- 
tions and data is fortuitous, since our calculations 
do not reproduce the shape of the ground state de- 
duced from the data. 

• For Po isotopes, the calculations overestimate the 
excitation energy of the 0^ level but reproduce the 
slope of the variation of the energy with N. For 
all isotopes for which data exist, our calculations 
predict the 0^ state to be oblate. 



2. Excitation energies 

The evolution of the energy of the first excited + 
states is compared to the experimental data in Fig. 12. 
The main findings are as follows: 

• Going from TV = 120 down to TV = 108, the 0^ and 
0jj~ excitation energies decrease rapidly in all the 
four isotopic chains, in agreement with the data. 

• As discussed above, the states 0^" and 0^ in Pb 
isotopes all have a well defined shape, the 0^" level 
being oblate above N = 106 and prolate below. 
The oblate and prolate configurations cross at N = 
106; in this nucleus the excitation energies of the 
O2 and 0j]" states are quite close and below 1 MeV. 

• The evolution of the energy of the 0^" state in 
Hg isotopes is in agreement with the data. In 
our calculations, this state can be called prolate 
only for N = 96, 98, 106 and 108; for all other 



3. Variation of mean-square charge radii 

The variation of the calculated mean-square charge 
radii is compared in Fig. 13 with the values deduced from 
the measured isotopic shifts. Our results for the Pb iso- 
topes obtained with the SLy6 and SLy4 Skyrme para- 
metrizations, each used with a different pairing strength, 
have already been compared with the data in Ref. [32]. 
In Ref. [18], results obtained with the SLy4 parametri- 
zation and two different values of the pairing interaction 
strength have been confronted with data. In both cases, 
it was observed that charge radii are extremely sensitive 
to the amount of mixing of mean-field wave functions 
corresponding to different shapes, which in turn is very 
sensitive to the pairing strength. 

A sudden increase of the charge radii is obtained at 
N = 104 in our calculation for the Hg isotopes. It is 
related to the shape transition seen in Fig. 3 and 9 that 
is, however, not corroborated by the data. By contrast, 
a jump of similar size has been observed in odd- A Hg 
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FIG. 13: (Upper panel) Variation of the charge radii <5(r 2 h ) for 
the lowest + states in Hg (normalized to the ground state 
(g.s.) of 194 Hg), Pb (normalized to the g.s. of 194 Pb), Po 
(normalized to the g.s. of 210 Po) and Rn (normalized to the 
g.s. of 204 Rn) isotopes, compared to the the experimental data 
for ground states taken from Refs. [18, 48]. 



[11, 13]. This observation indicates that two states based 
on either prolate or oblate shapes indeed coexist but that 
our calculations do not predict correctly their relative po- 
sition. Above TV = 104, the ground-state wave functions 
are weakly oblate, but become strongly peaked at prolate 
shapes for 100 < TV < 104, as shown in Fig. 9. However, 
this onset of large prolate deformation in the ground state 
is in contradiction with the experimental data. Clearly, 
the dominant configuration of the ground state should 
not change much down to at least TV = 102. 

In the Rn isotopes, the kink at TV = 112 in the calcu- 
lated values of 6(r 2 h ) is due to the onset of weakly oblate 
deformations in the ground state at TV — 110, cf. Figs. 9 
and 12. 



4- Monopole transition strength 

Figure 14 displays the E0 strength p E0 , Eq. (10), be- 
tween the 4 " and the 4 " states. This electric monopole 
transition strength is correlated to the size of the defor- 
mation and to the amount of mixing between configura- 
tions corresponding to different shapes [2]. In general, 
large EQ strength arises from a strong mixing between 
states with different radii. 

Compared to the very sparse experimental data that 
exist for some Hg, Pb and Po isotopes, our calculated 
p E0 values are too large by an order of magnitude. Ac- 
cording to Rcf. [2], this indicates a too stong mixing of 
configurations with different deformation in at least one 
of the two 4 " states. 
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FIG. 14: (Color online) Electric monopole transition 
strengths p% for the transitions between the first excited + 
state and the ground state (cf. Fig. 12) in Hg, Pb, Po and 
Rn isotopes. Experimental data are taken from Ref. [49]. 



F. Low-lying states in the TV = 106 isotonic chain 

After the discussion of the J = states, we will now 
turn to the properties of states with finite angular mo- 
mentum. Before looking at their systematics, let us first 
analyze how these levels group into bands in the case 
of the TV = 106 isotones. For this neutron number, our 
calculations predict three nearly-degenerate + levels for 
Hg, Pb and Po, cf. Figs. 8-10, which makes these isotopes 
the ideal laboratory for this purpose. 

The projected energy curves arc plotted in Fig. 15 for 
186 Hg, 188 Pb, and 190 Po. The energy of the correlated 
GCM states are indicated by bars located at their mean 
deformation /?,/&. The energy curves projected only on 
particle numbers are also provided for comparison. There 
is no rotational band on top of the ground state of 188 Pb, 
which confirms our interpretation that this level is a cor- 
related spherical state. However, the mean deformations 
(3jk and energies of the GCM states, as well as the sys- 
tematics of B{E2) values and spectroscopic quadrupole 
moments to be discussed below, indicate two rotational 
bands in this nucleus, one of prolate and the other of 
oblate deformation. Similar bands are found for 186 Hg 
and 190 Po. For both bands and in all three isotones, the 
+ band head is strongly perturbed. 

The energy of the + band head is even pushed above 
that of the 2 + state for the oblate band in Hg and Po and 
quasi-degenerate in Pb. For larger- J values, the bands, 
however, become more regular with a steady increase of 
the mean deformation j3jk as a function of J. This means 
that deformation follows closely the minimum of the pro- 
jected energy curve for the oblate bands, whereas it is 
becoming larger for the prolate sequences. In any event, 
a third low-lying + state that cannot be associated with 
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FIG. 15: (Color online) Deformation energy curves projected 
onto particle number (N&Z), or additional onto angular mo- 
mentum (J = 0, 2, . . . , 8), as well as the projected GCM 
states in the N = 106 186 Hg, 188 Pb, and 190 Po isotones. The 
projected energy curves are plotted as a function of the in- 
trinsic deformation ft of the mean field states. The energies of 
projected GCM states are indicated by bullets and horizontal 
bars placed at the average deformation pjk- All energies are 
normalized to the + ground state. 



a rotational band is also found for 186 Hg and 190 Po. 

The theoretical spectra of the three nuclei are com- 
pared with the experimental data in Fig. 16. As al- 
ready discussed in previous studies [31], the calculated 
spectra are too much expanded. However, the overall 
trends of experiment are rather nicely reproduced. The 
intraband E2 transition strength B(E2 : J + 2 — > J) in- 
creases steadily with angular momentum, which is consis- 
tent with the picture shown in Fig. 15 that the collective 
wave functions are mainly built from either purely pro- 
late or oblate mean-field states as soon as J > 2. Both 
the experimental B(E2) values within and out of bands 
are reproduced satisfactorily by our calculation. 



G. Systematics of 2 + and 4 + states 

The variation as a function of TV of the excitation ener- 
gies of the low- lying + , 2 + and 4 + states are presented in 
Figs. 17-20 for the Hg, Pb, Po and Rn isotopes. The ex- 
perimental data are taken from Rcfs. [16, 46, 47, 50, 51]. 
Lines join levels with a similar structure in two successive 
isotopes. 

As already discussed for the N = 106 isotones, the 
spectra are expanded too much. There can be several 
sources to this deficiency: a variational space that is bet- 
ter suited to optimize the + ground state than the ex- 
cited states, triaxiality effects, the treatment of pairing 
correlations or deficiencies of the energy density func- 
tional. We will focus on the evolution of the relative 
position of the rotational bands as a function of N (and 
sometimes of Z), for which the spacing between levels in 
a band is in general not crucial. The main findings are 
as follows: 

• The experimental trends are reasonably well re- 
produced by our calculations for the Hg and the 
Pb isotopic chains. In Pb, the prolate and the 
oblate bands cross between N — 108 and N = 106, 
in agreement with the data. For the Hg iso- 
topes, there is an exchange of structure between 
the ground state and the first excited + level in 
such a way that we obtain a prolate ground state 
at N = 104 to 100, in contradiction with the data. 
However, the first two + states remain close in 
energy for these three isotopes. 

• In Hg and Pb, the mixing between the + levels 
is usually too large in the calculations, distorting 
too much the rotational character of bands at low 
spin. Experimentally, the energy of the first 2 + 
level is nearly constant over the entire isotopic Hg 
chain, whereas its variation is much larger in our 
calculations. The energy difference between the 2 + 
and 4 + levels varies in a manner more similar to 
the data than the excitation energy of the 2 + level. 

• In the Hg and Pb isotopic chains, the calculated 
excitation energies of the 2f and 4^~ levels reach 
their lowest value at N — 102. This result agrees 
with the experimental data for the Hg isotopes. For 
the Pb isotopes, the lowest excitation energy for 
these levels is obtained in the experimental data at 
N = 104, but the calculated results at N = 102 
and 104 are only marginally different. 

• There is a clear disagreement between calculations 
and experimental data for the heavy Po isotopes. 
The overestimation of energies is close to a factor 
of two and much larger than for other nuclei. In 
particular, the near-constant value of the experi- 
mental 2 + level energy between N = 120 and 112 
is not reproduced at all. Being close to the N = 126 
shell closure, these isotopes are probably not collec- 
tive and would require to include np-mh excitations 
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FIG. 16: (Color online) Comparison between the calculated and the measured low-lying spectra for 186 Hg, 188 Pb, and 190 Po. 
The B(E2) values are given in e 2 b 2 units. The experimental data are taken from Refs. [8, 17, 52-54]. 
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FIG. 17: (Color online) Systematics of calculated (upper 

panel) and experimental (lower panel) low-lying + ,2 + and FIG. 18: (Color online) Same as Fig. 17, but for the Pb iso- 
4+ states in the Hg isotopes. The experimental data for the topes, 
lowest 6 + and 8 + states are given as well to indicate the evo- 
lution of the band structure. 
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FIG. 19: (Color online) Same as Fig. 17, but for the Po iso- 
topes. 
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FIG. 20: (Color online) Same as Fig. 17, but for the Rn iso- 
topes. 
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FIG. 21: The quadrupole deformations /3 (s) (cf. Eq.(8)) of 
the two lowest J = 2 states calculated from the corresponding 
spectroscopic quadrupole moments as functions of the neutron 
number. 



that are not generated by a constraint on deforma- 
tion. The calculated spectra are more realistic for 
the isotopes with 116 < N < 106. 

• The disagreement with the data is not as strong 
for the Rn isotopes. In particular, the energy dif- 
ference between the 2+ and 4 + levels is closer to 
the experimental value. 

The quadrupole deformation parameters determined 
from the spectroscopic quadrupole moment, (3^ (cf. 
Eq.(8)), are plotted as a function of the neutron number 
in Figs. 21 and 22 for the lowest two 2+ and 4 + states, 
respectively. This parameter, which is directly related to 
an observable, permits a much less ambiguous assignment 
of a prolate or an oblate character to the levels than the 
mean deformation /3jk of the wave functions. A crossing 
between oblate and prolate configurations takes place at 
N = 106 for the Hg, Pb and Po isotopes, at N = 108 for 
Rn. Again, this result is not confirmed by the data for 
the Hg isotopes. For the heaviest Po and Rn isotopes, 
/?( s ) is very small. This is a sign that the 2 + and 4 + 
are probably not based on a deformed state, but would 
be better described by np-mh excitations of the kind al- 
ready mentioned above. This could also be at the origin 
of the much too large excitation energies that we have 
obtained for these 2 + and 4 + states (see Fig. 19 and 20); 
our purely collective GCM basis not being well-suited to 
describe these states. 

The quadrupole deformations (3^(2^ — > 0i), cf. 
Eq.(9), corresponding to the electric quadrupole transi- 
tion strengths B(E2 : 2% — > Ojfc) are plotted in Figure 23 
for the transitions from the two first 2 + to the + ground 
state and the Oj level, together with the available experi- 
mental data. The assignment of oblate and prolate labels 
to the 2 + states are based on spectroscopic moments and 
are, thus, non ambiguous. The main findings are as fol- 
lows: 
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FIG. 23: (Color online) Systematics of the deformations 
derived from the reduced E2 strengths for the transition from 
the predominant prolate and oblate 2 + states to the first two 
+ states as a function of neutron number. The data for 
transition from the 2f state to the ground state are taken 
from Refs. [8, 14, 55-58]. 



• The systematics of (3^ values for transitions be- 
tween oblate or prolate states are quite similar in 
the four isotopic chains. 

• For i80.i82,i84 Hg; the experimental /3 (t) (2+ -> 0+) 
values are in good agreement with the calculated 
ones for the transition from the oblate 2 + level to 
the oblate 0^" state. It shows also that theoretical 
transition probabilities can be accurate in spite of 
the fact that the energies and the order of levels are 
not reproduced. 

• For the Pb isotopes with N = 104 and 106, the 
experimental (3^(2^ — > Of) value is in good agree- 
ment with the calculated transition from the oblate 



2 + state to the spherical 0^ state. It suggests 
that, with decreasing neutron number, the 2^ state 
keeps predominantly an oblate character down to 
N = 104, beyond which the prolate 2 + state be- 
comes the yrast state in the calculations. 

• For the N = 116 Po isotope, the experimental 
(3^(2^ — > 0+) value is reproduced quite well by 
the calculation. It corresponds to a transition 
from the oblate 2f level to the spherical Of state. 
As the neutron number decreases to N = 112 and 
110, the transitions calculated to the oblate 0^" 
state are in better agreement with the experimen- 
tal /3^(2^ —> 0^) values than the transitions to 
the spherical Of state. This indicates an onset of 
oblate deformation already at N = 112, which is 
consistent with the observation made from the ex- 
perimental isotope shifts, c.f. Fig. 13. 

• There are no experimental data for the B(E2) val- 
ues in the Rn isotopes discussed here. 



H. Kinetic moment of inertia along the yrast band 

The variation of the moment of inertia along a band 
can be an indicator of its nature. There are well estab- 
lished generic properties of the moments of inertia such 
as their increase with deformation, or that for the same 
value of |/3 1 the moment of inertia is larger for prolate 
deformations than for oblate ones. The kinetic moment 
of inertia Q^ 1 ) is defined as 



Cj(D 



(J) = 



hy/j{J+l) 

w(J-> J-2) 



and the frequency w, 
hw(J -> J - 2) = 



[E x {J)-E x {J-2)\ 



(13) 



(14) 



The kinematic moments of inertia for transitions along 
the yrast line are plotted in Fig. 24-27 for transitions 
from J = 2, 4, . . . , 10 to J - 2 in Hg, Pb, Po and Rn 
isotopes, respectively. As discussed above, the calcula- 
tions underestimate the moment of inertia, cf. Fig. 16. 
To make the comparison with the experimental data eas- 
ier, we have divided the calculated rotational frequency 
by a factor 1.50. 

There are several common features in the four isotopic 
chains: 

• For a static rotor, the kinematic moment of inertia 
is constant. The only nuclei in our sample which 
came close to this idealized case are the Hg isotopes 
with N — 112 and 114, at least up to J = 8 in both 
calculation and experiment. For all other nuclei the 
pattern is more complex. 
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FIG. 27: Same as Fig 24, but for the Rn isotopes. 



• The kinematic moment of inertia for J = 2 
is significantly different from those corresponding 
to transitions from higher spins. This can be at- 
tributed to the ground state configuration that in 
most cases is dominated neither by an oblate nor 
by a prolate configuration. In the Pb isotopes, the 
ground state is a correlated spherical state, with a 
very different structure from the first 2 + state. In 
any case, the ground state cannot be considered as 
the band head of a rotational band. 

• The kinematic moments of inertia 3^ of states 
with spin J > 4 can be divided into two groups: 
the isotopes with N < 106 have a prolate configu- 
ration and a larger SyW, whereas the isotopes with 
N > 108 arc oblate and have a smaller SjW. This 
feature is present in both the calculations and the 
data. This confirms again that prolate configura- 
tions dominate the yrast states below N = 108. 

• For many nuclei, the calculated Q^ 1 -* moments take 
almost constant values for the highest angular mo- 
menta studied here. Such behavior, however, is al- 
most never seen in the data. This discrepancy can 
be expected to have the same origins as the overall 
expansion of the spectra, cf. Sect. IIIG. 



IV. SUMMARY AND OUTLOOK 

We have shown that a detailed comparison between 
the results obtained by a beyond mean-field method and 
experimental spectroscopic data can be performed for a 
large number of nuclei. This comparison can be carried 
out without ambiguities when restoring quantum num- 
bers and thereby selection rules for transitions. The re- 
gion around the neutron-deficient Pb isotopes is of par- 
ticular interest for such studies because it presents a very 
rich variety of phenomena, and the structure of these nu- 
clei varies rapidly as a function of the neutron number. 
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The results of our model are in qualitative agreement 
with the experimental data. We obtain the coexistence 
of several low-lying + levels showing mixing between 
oblate, prolate and spherical configurations. The order 
of these levels is sometimes different from that deduced 
from the analysis of the data for a-decay hindrance, radii, 
moments of inertia and electromagnetic transition den- 
sities. However, the theoretical B(E2) values within a 
rotational band of given prolate or oblate nature are usu- 
ally described well even though the relative position of 
the band heads is not correct. 

We find that labeling states as oblate or prolate might 
not always be adequate, especially for nuclei in this mass 
region. In particular, the + and 2 + states often re- 
sult from a mixing between symmetry-restored mean- 
field wave functions corresponding to shapes spread over 
a wide range of deformations. The role of the mixing di- 
minishes with increasing angular momentum, such that 
the spectra have more apparent rotational or vibrational 
character at high spin. 

Our study has been limited to the mixing of particle- 
number and angular-momentum restored time-reversal 
invariant axial configurations. Tools to extend such cal- 
culations to triaxial shapes have been set-up recently [26- 
28, 59-61]. However, the presently available computa- 
tional resources do not yet allow for their systematical ap- 
plication to very heavy nuclei. We have seen that triaxial 
configurations might play a role for some of the lightest 
Hg and some of the heaviest Po and Rn isotopes studied 
here. The absence of this degree of freedom is one of the 
reasons for the systematically too expanded theoretical 
excitation spectra. To further improve the description 
of excited states, a second extension of the variational 
space to be considered in the future is the projection 
and mixing of time-reversal invariance-breaking mean- 
field states. On the one hand, starting from cranked 
HFB states that are optimized for finite angular mo- 
menta will incorporate the rotational alignment of single- 
particle states and the weakening of pairing correlations 
at high spin that comes with it. And on the other hand, 
the use of so-called blocked HFB states of the broken-pair 
type will provide a natural starting point to describe the 
non-collective states that can be suspected to dominate 
the low-lying excitation spectra of the heaviest Po and 
Rn isotopes discussed here. Developments in both direc- 
tions are underway, but, again, first applications will be 
limited to light nuclei for computational reasons. 

The energy density functionals currently used do not 
permit yet to reach an accurate quantitative description 
of spectroscopic data. However, we recall that the en- 
ergy density functional used here has quite a simple form 



and that its parameters have been adjusted on the mean- 
field level to bulk properties of nuclear matter and magic 
nuclei. Hence, the overall good qualitative description of 
the rapidly evolving collective states of the nuclei studied 
here should be viewed as a success. This can be partially 
attributed to the generic features of the single-particle 
level schemes that remain an efficient tool to qualitatively 
analyze also results from bcyond-mean-ficld calculations, 
as has been shown by our discussion. In addition to the 
well-known correlation between the opening and closing 
of gaps around the Fermi energy in the Nilsson diagram 
and the appearance of the various minima in the deforma- 
tion energy surfaces, we have also seen that the spread of 
the correlated collective wave functions over different de- 
formations is correlated with the deformations at which 
the intruder levels cross the single-particle levels occu- 
pied in the spherical configuration. In general, we find 
that the mixing of states with a different number of oc- 
cupied intruder levels is disfavored, which is consistent 
with the usual classification of shape coexisting states in 
the interacting shell model. A better quantitative de- 
scription of fine details of shape coexistence, such as the 
relative energy between oblate and prolate band heads 
in the Hg isotopes below N = 108, will require improved 
single-particle energies, cf. also the discussion in Ref. [62]. 
First explorative studies in this direction indicate, how- 
ever, that it is unlikely to obtain a significantly better 
description of single-particle energies by a refit within the 
current standard form of the energy density functionals 
[63, 64], and that extensions are needed. Work in this 
direction is also underway. 
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